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ABSTRACT 

We present the combined results of three photometric monitoring campaigns targeting 
very low mass (VLM) stars and brown dwarfs in the young open cluster IC4665 (age 
^ 40 Myr) . Each of our observing runs covers timescales of ~ 5 days in the seasons 
1999, 2001, 2002, respectively. In all three runs, we observe ~ 100 cluster members, al- 
lowing us for the first time to put limits on the evolution of spots and magnetic activity 
in fully convective objects on timescales of a few years. For 20 objects covering masses 
from 0.05 to 0.5 Mq we detect a periodic flux modulation, indicating the presence of 
magnetic spots co-rotating with the objects. The detection rate of photometric periods 
(~ 20%) is significantly lower than in solar-mass stars at the same age, which points 
to a mass dependence in the spot properties. With two exceptions, none of the objects 
exhibit variability and thus spot activity in more than one season. This is contrary 
to what is seen in solar-mass stars and indicates that spot configurations capable of 
producing photometric modulations occur relatively rarely and are transient in VLM 
objects. The rotation periods derived in this paper range from 3 to 30 h, arguing for a 
lack of slow rotators among VLM objects. The periods fit into a rotational evolution 
scenario with pre-main sequence contraction and moderate (40-50%) angular momen- 
tum losses due to wind braking. By combining our findings with literature results, we 
identify two regimes of rotational and magnetic properties, called C- and I-sequence. 
Main properties on the C-sequence are fast rotation, weak wind braking. Ha emission, 
and saturated activity levels, while the I-sequence is characterised by slow rotation, 
strong wind braking, no Ha emission, and linear activity-rotation relationship. Rota- 
tion rate and stellar mass are the primary parameters that determine in which regime 
an object is found. We outline a general scheme to understand rotational evolution 
for low-mass objects in the context of these two regimes and discuss the potential as 
well as the problems of this scheme. 

Key words: stars: low-mass, brown dwarfs, stars: rotation, stars: evolution, stars: 
activity 



1 INTRODUCTION 

Rotation and magnetic activity are key parameters of stellar 
evolution. While these parameters have been well-explored 
for solar-mass stars over the past decades, our knowledge is 
still sparse in the very low mass (VLM) regime. The inter- 
est in rotation and activity in VLM stars and brown dwarfs 
is fueled by the hypothesis that these properties depend on 
a fundamental level on the interior structure. Specifically, 
to operate the type of dynamo that generates the solar- 
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type cyclical magnetic field, the presence of a shear layer 
between convective envelope and ra diative core is thought 
to be of prime importance (e .g. ISpiegel fc Weiss! Il980l : 
ISpruit fc van BallegooiienI 1 19821 . and references therein). 
Objects with masses < 0.3 Mq, however, are fully convec - 
tive throughout their evolution (|Chabrier fc Baraffelll997l ). 
Thus, if the aforementioned assumptions are correct, they 
would not be able to generate magnetic fields in the same 
way as solar-type stars. Changes in magnetic properties and, 
consequently, rotational evolution are therefore expecte d at 
the fully-convective boundary (see iDelfosse et a ll ll998h . 

Rotation and activity are interdependent in a twofold 
sense: On one side, rotation plays an important role in driv- 
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ing the dynamo and essentially powers magnetic activity. 
On the other side, stellar winds driven by magnetic fields 
are the predominant agent of the spindown on the main- 
sequence and thus determine the long-term rotational evo- 
lution of stars. For solar-mass stars, the close connection 
between rotation and activity manifests itself in the rota- 
tion/activity relation, which is linear for slow rotators and 
flatt ens at high rotat ion rates, as well as in the Skumanich 
law (|Skumanichll 197^ ) . which establishes that rotation rates 
and activity indicators both decline with the square root of 
age on the main-sequence. 

Over the past decade, it has been shown convincingly 
that, with the exception of a few very young objects, VLM 
stars and brown dwarfs are fast rotators with periods 
gener all y shorter than 2 d (e.g. iBailer- Jones fc Mundtl 



2004bl: 'Herbst et al 



2OOII: IScholz fc Eislofel_ 

Reiners fc Basril ~ 20081 : iRodrfguez-Lc dcsma et al. 
The periods found in the Monitor project confirm this 
finding (see in particular llrwin et all bOOSbT ). Rotational 
braking due to stellar winds is clearly less efficient in this 
mass regime. On the other hand, a number of ultracool 
objects are now known to harbour strong, large-scale 
magnetic fields, as evidenced by Zeeman Doppler imaging 
l|Donati et akl l2008l: iMorin et al.l 12008^ ). direct field mea- 



sure ments ( Reiners fc Basri 
(e.g. lBergerll2006l : 



2007 



and radio observations 

Hallinan et"all 120061 ). A full theory that 



accounts for all these findings is still missing. 

Investigating rotation and activity in the VLM regime 
therefore has the potential to probe fundamental physical 
questions: To what extent does interior structure affect the 
observable properties of stars? How does magnetic field gen- 
eration change as a function of object mass? Can we suc- 
cessfully predict rotation rates for a given age and mass? 
The latter problem leads to the intriguing possibility to use 
rotation periods - a quantity that can be measured with 
high accuracy - to determine stellar ages and masses (gy- 
ro chronology, see iBarnedlioOTl ). 

The goal of this paper is to provide new observational 
constraints on rotation and activity in the very low mass 
regime. We will probe a specific subset of magnetic phenom- 
ena, namely the properties of magnetically induced spots in 
the photosphere. We present results from three photometric 
monitoring campaigns in the young open cluster IC4665, in- 
cluding 20 new rotation periods for VLM stars and brown 
dwarfs (Sect. I2I4I1 . Through repeated monitoring, we are 
able to probe the characteristics and the long-term evolution 
of magnetic spots in this critical mass regime (Sect. [5]). Fi- 
nally, we will use our periods in combination with literature 
data to probe rotational braking and magnetic properties in 
low-mass objects (Sect. [6] and [7]). 



2 PHOTOMETRIC MONITORING AND DATA 
REDUCTION 

2.1 Observations 

The target field for this variability study covers parts of the 
central region of the open cluster IC4665. The field is centred 
at 17''45'"15.0^5 05°21'17'.'0 (J2000.0) and is covered by 
our photometric survey (Eisloffel et al., in prep., see Sect.|3]). 
This region was observed in three I-band monitoring cam- 
paigns covering a time span of four years in total, using the 



Table 1. Observing log for the three observing runs, runs A and C 
were carried out with the ESO/MPG 2.2-m WFI, run B with the 
TLS Schmidt camera. Columns 3-6 contain number of exposures, 
integration time per exposure, sky conditions, median seeing. 





date 




(Ci Yl~l tim*::* 






A 


28/05/1999 


1 


600 s 


photometric 


1"2 




05/06/1999 


44 


500 s 


some cirrus 


0'.'7 




06/06/1999 


43 


500 s 


photometric 


o'.'g 




08/06/1999 


3 


500 s 


partly cloudy 


1'.'2 




09/06/1999 


5 


500 s 


photometric 


1'.'2 


B 


19/05/2001 


13 


600 s 


partly cloudy 


3'.'4 




20/05/2001 


24 


600 s 


some cirrus 


2" 7 




21/05/2001 


23 


600 s 


partly cloudy 


2'.' 9 




22/05/2001 


23 


600 s 


photometric 


2'.' 6 




23/05/2001 


22 


600 s 


photometric 


2'.' 6 


C 


29/05/2002 


15 


500 s 


photometric 


1'.'6 




07/06/2002 


14 


500 s 


partly cloudy 


1'.'2 




08/06/2002 


14 


500 s 


photometric 


I'.'S 




09/06/2002 


10 


500 s 


some cirrus 


1'.'5 




12/06/2002 


17 


500 s 


some cirrus 


1'.'6 




13/06/2002 


17 


500 s 


some cirrus 


I'.'S 



ESO/MPG 2.2-m Wide Field Imager at La Silla and the 
2-m Schmidt telescope at the Thiiringer Landessternwarte 
(TLS). Both instruments have a comparable field of view of 
~ 0.32 sqdeg. 

Details for the three observing runs are given in Ta- 
ble [T] In the first WFI run in May/June 1999, we took 96 
long exposure images of this field, including complete cov- 
erage of two consecutive nights. The second WFI run was 
carried out in service mode, resulting in dense coverage of 
typically 2-3 h per night for six nights distributed over two 
weeks. This campaign provided 87 time series images of our 
field in IC4665. The TLS observing run covers five consecu- 
tive nights with typically 4 h of continuous observations per 
night, in total 105 exposures. The TLS time series suffered 
from bright sky background and sub-optimal seeing, and 
thus the photometric accuracy was significantly lower than 
in the WFI runs. The typical integration time for individual 
exposures was 500 sec in all three runs. 



2.2 Image reduction 

All images were re duced in a similar fashion, follo wing the 
recipes outlined in lScholz fc Eisloffell (|2004al . 120051 ). The re- 
duction includes bias subtraction and flatfield correction us- 
ing either twilight fiats (WFI) or domeflats (TLS). Addi- 
tionally, the I-band images have to be corrected for multi- 
plicative large-scale illumination gradients, a consequence 
of imperfect flatfielding, and additive, small-scale fringe 
structures. Both effects can be corrected using a 'superflat' 
constructed from images of dark sky regions. By median- 
filtering the superfiat, large-scale structures (illumination) 
can be separated from the small-scale fringes. As a result, we 
obtain an image of the average illumiation and the average 
fringe pattern. Since fringes are variable as a function of air- 
mass, temperature, and weather conditions, an appropriate 
scaling is required for each time series i mage. We used the 
autom ated scaling procedure outlined in IScholz fc Eisloffell 
(|2005l ). and manually refined the resulting scaling factors if 
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Figure 1. RMS vs. relative magnitude for all three runs: solid lines show the binned median rms, large symbols the VLM candidate 
members. Note the different scale on the y-axis of the middle panel. For run B, we only plot 4000 of 10000 field objects, for clarity. Due to 
the differences in the I-band filters used in the various campaigns, the relative magnitudes are not easily transformed to calibrated I-band 
photometry. The approximate I-band mags for the zeropoint on the relative scale are 17.2, 16.2, and 16.5 mag for runs A-C, respectively, 
with an uncertainty of ±0.3 mag. These numbers are derived by comparing the relative magnitudes with the absolute photometry derived 
in our survey (Sect. [Sjl. 



necessary. We estimate the residuals after fringe correction 
to be not larger than 1% of the sky background. 



plotted in large symbols in Fig. [T] are the datapoints for the 
sample of cluster member candidates, described in Sect.|31 



2.3 Instrumental magnitudes 

We derived instrumental magnitudes for all objects in the 
time series field by fitting their point spr ead function ( PSF) 
using the daophot package within IRAF l|Stetsonlll987h . For 
the WFI runs, all photometry steps were done on each chip 
of the mosaic separately to account for the variable proper- 
ties of the CCDs, e.g. the nonlinearity of the chip response. 

Model PSFs were calculated using typically 40 isolated 
stars. Then, this model was fit to each object in the field. 
The reliability of the resulting photometry was checked by 
subtracting all PSF fits from the original image. As result we 
obtained essentially empty frames, in which only saturated 
objects, bright galaxies, and some close undetected binary 
companions remained. 



2.4 Relative calibration 

The lightcurves of all objects were calibrated relatively with 
an average lightcurve from a set of non-variable reference 
stars in the observed field. The selection of reference stars 
was performed in a similar fashi on for all three campaigns , 
using the procedure described in IScholz fc Eisloffell (l2004al ). 
In short, the routine compares each potential reference star 
with the average of all other potential reference stars. For 
the WFI runs, this procedure was again carried out for each 
chip of the mosaic separately. Typically, the final sample of 
reference stars comprises of ~ 100 objects per chip. From 
these objects, the final reference lightcurve was determined. 
By subtracting this reference lightcurve from all time series 
from the respective observing run, we obtained relatively 
calibrated magnitudes for all objects. 

In Fig. [1] a-c we plot the lightcurve RMS vs. relative 
magnitude for all three runs. Overplotted is the median rms 
per magnitude bin as solid line after excluding outliers - this 
function gives the typical photometric error at a given rel- 
ative magnitude. For bright, unsaturated stars with I-band 
magnitudes of 15-16, we obtain average uncertainties of 5- 
8mmag in runs A and C and 10-20 mmag in run B. Ovcr- 



3 THE SAMPLE OF PROBABLE CLUSTER 
MEMBERS 

The targets for our variability analysis are 113 sources clas- 
sified as probable cluster members in our deep photomet- 
ric survey in IC4665. The detailed outcomes of this survey 
will be published in a forthcoming paper (Eisloffel et al., 
in prep.). In short: We surveyed ~ 0.4 sqdeg of the central 
cluster in the R- and I-band using the ESO/MPG WFI at 
the 2.2-m telescope at La Silla - these images have been ob- 
tained in the same observing run that provided the first time 
series (run A, see Sect. I2.l"|) . The R- and I-band data is cali- 
brated by comparing with photome try of Landolt standard 
fields observed in the same nights (lLandoltlll992l ). A deep 
(I,R-I) colour magnitude diagram yielded a sample of 118 
objects whose optical colours are consistent with an empir- 
ically determined isochrone for the cluster IC4665 (see Fig. 
[2}. We complemented the survey with near- infrared pho- 
tometry, either from 2MASS or from our own J-band dat4i]- 
From the primary sample, 112 objects show near-infrared 
colours in ag reement with t he mo del isochrone for an age of 
40Myr from lBaraffe et all (119981 ). 

Furthermore, we surveyed a larger area of the clus- 
ter with the 2-m Schmidt telescope at the Thiiringer Lan- 
dessternwarte Tautenburg in a run in May 1999. The I-band 
photometry from this campaign combined with 2MASS 
near-infrared data yielded an (I,I-J) colour-magnitude di- 
agram. Applying a colour selection, we obtained an addi- 
tional sample of 57 probable cluster members. Only one of 
these sources is located in the monitoring field and shows 
a significant period in run A of the variability study (ob- 
ject #119, see below). We include this object in the discus- 
sion. Thus, the sample for the time series analysis comprises 
113 objects, with calibrated photometry in optical and near- 
infrared bands. 



^ Observations carried out with the near-infrared camera Q' at 
the 3.5-m telescope on Calar Alto/Spain, Deutsch-Spanisch As- 
tronomisches Zentrum 
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Table 2. Basic data (coordinates, photometry, masses) for the cluster member candidates with periodic lightcurve in at least one of 
the monitoring campaigns. Full survey results to be published in Eisloffel et al. (in prep.). "2M" indices designate photometry from th e 
2MASS database. Mass estimates derived by comparing the near-infrared photometry with evolutionary tracks bv lBaraffe et al 
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09.7 


17.59 


15.74 


13.90 


13.93 


13.31 


13.00 


0.51 


40 


17 


44 


46.45 


5 


32 


01.9 


20.04 


17.83 


16.16 


15.77 


15.28 


14.89 


0.14 


49 


17 


44 


54.68 


5 


23 


52.4 


17.71 


16.07 


14.24 


14.49 


13.78 


13.56 


0.37 


55 


17 


45 


05.57 


5 


11 


11.0 


18.11 


16.34 


14.67 


14.67 


14.01 


13.78 


0.32 


57 


17 


45 


06.69 


5 


30 


38.4 


18.69 


16.79 


15.28 


15.07 


14.41 


14.01 


0.25 


59 


17 


45 


10.95 


5 


33 


28.5 


19.14 


17.14 


15.53 


15.43 


14.67 


14.49 


0.19 


63 


17 


45 


13.53 


5 


18 


39.2 


17.90 


16.08 


14.33 


14.32 


13.73 


13.36 


0.41 


73 


17 


45 


21.39 


5 


34 


37.8 


19.39 


17.32 


15.73 


15.48 


14.80 


14.59 


0.18 


75 


17 


45 


26.00 


5 


29 


40.5 


20.98 


18.48 


16.56 


16.21 


15.63 


15.24 


0.11 


76 


17 


45 


27.20 


5 


33 


57.8 


23.00 


20.32 


18.00 








0.05 


88 


17 


45 


37.44 


5 


23 


40.5 


19.44 


17.63 


15.88 


15.95 


15.27 


15.00 


0.13 


100 


17 


45 


49.89 


5 


09 


40.2 


18.23 


16.32 


14.61 


14.66 


13.97 


13.64 


0.34 


110 


17 


45 


59.35 


5 


29 


31.6 


20.57 


18.24 


16.32 


16.01 


15.44 


15.12 


0.12 


118 


17 


46 


06.29 


5 


36 


39.3 


20.87 


18.50 










0.09 


119 


17 


45 


37.10 


5 


23 


08.1 


15.79 






14.28 


13.63 


13.34 


0.42 




0.5 1 1.5 2 2.5 

R-I (mag) 



Figure 2. I, R-I colour magnitude diagram for our time series 
field in IC4665: The two solid lines illustrate the colour criterium 
applied to select the primary sample. All probable cluster mem- 
bers according to their colours are shown as large dots; objects 
with periodic lightcurve are marked with squares. Errorbars indi- 
cate typical photometric uncertainties. For clarity, only ^ 10000 
out of ~ 30000 total datapoints are shown. Photometry obtained 
during run A (see Table [Til. Missing from that figure is object no. 
119, which was identified in a different campaign (see text). 



The photometrically selected cluster member candi- 
dates may be contaminated by red field objects. The con- 
tamination rate was estimated based on optical colour- 
magnitu de diagrams sim ulated usin g the Galaxy model pre- 
sented bv lRobin fc Creze (WSG) and lRobin et all l|2003t ). We 
derive a contamination rate of ~ 30% for the total candidate 
sample, most of the contaminating objects are expected to 
be M dwarfs in the fore- or background of the cluster . From 
the space density for M dwarfs given bv lReid fc Gk3 (|l997l ) 



and our survey area follows a contamination rate of about 
30%, consistent with the Galaxy model. 

In Fig. [2] we show the (I, R-I) colour magnitude dia- 
gram, probably cluster members are marked. In our three 
monitoring campaigns, we found 20 objects with significant 
periodicity in the lightcurve (see Sect.|3|; these objects are 
marked with open circles. As can be seen in this plot, the 
majority of the periodic objects forms a tight sequence with 
less scatter than the total sample, as expected for objects on 
the cluster isochrone. This may be seen as evidence for lower 
contamination among the objects with period, although for 
a definite decision spectroscopy is required. 

In Table [2] we provide the basic data for all objects 
in our survey with period detection in at least one of the 
monitoring campaigns. We estimate masses for these ob- 
jects by comparing their near-infrared magnitudes with the 
cluster isochrone from Baraffe et al. (1998), assuming an 
age of 40Myr, a distance of 350 pc, and an extinction of 
£s-v = 0.18 (|Mermilliodlll981al lbh. According to this esti- 
mate, the objects span a mass range from 0.05 to O.5Af0 and 
thus include a few brown dwarf candidates. Due to uncer- 
tainties in the theoretical evolutionary tracks and differences 
between the standard bands and the filters used in the sur- 
vey, resulting in calibration offsets, the mass estimates may 
be inaccurate by as much as 100%, particularly at very low 
masses. Since both effects are systematic, the masses for our 
objects are comparable on a relative scale. 



4 TIME SERIES ANALYSIS 

The main focus of the time series analysis was the search for 
periodic variabilities. The period search is described in detail 
in Sect. 14. Il In addition, we looked for conspicuous signs for 
clear non-periodic variability in the lightcurves: All candi- 
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date lightcurves were inspected visually to register obvious 
signs of rapid variability. In particular, we searched for signs 
of flares, i.e. for brightness eruptions with gradual decline 
lasting several hours. No such events with > 0.1 mag ampli- 
tude were found, indicating that flare events are rare in VLM 
objects. Based on our time coverage, we estimate an I-band 
flare rate of < 3.2-10"'* , in agreement with the result ob- 
taine d in a previous paper (2.3- 10-* /t~ MScholz fc Eisloffell 
l2005l ). 

In the RMS vs. magnitude plots in Fig. [T] we mark the 
cluster members with larger symbols, to allow an assessment 
of the generic variability. As can be seen in these plots, VLM 
objects in IC4665 show only low-level photometric varia- 
tions, with amplitudes < 0.05 mag. Most of the objects that 
are observed to have an RMS significantly larger than the 
field average are later identified as objects with periodic vari- 
ability. This is in contrast to the findings for field L dwarfs 
with similar masses, where variable objects are rarely found 
to exhibit a photometric period (e.g., iBailer- Jones fc Mundtl 
l200ll : lGelino et al.ll2002l ). 

4.1 Period search for runs A and C 

We used the Scargle periodogram l|Scarglelll98'3 ) to identify 
candidate periods in the lightcurves from runs A and C. For 
the initial detection, we relied on a preli minary significance 
estimate based on the equation given bv lHorne fc Baliunad 
l|l986t ). All candidate periods are then verified in a num- 
ber of independent tests. In particular, periodograms and 
lightcurves of possible periodic objects were compared with 
those of nearby stars and excluded if the neighbours show 
a similar periodicity. In addition, we CLEANed the peri- 
odograms for all periodic objects from sam pling artefacts 
and a liases using the procedure given by 'Robe rts et al.l 
([l983) 

and checked if the peak found in the Scargle peri- 
odogram is still present. Only periods which pass all tests are 
ac cepted. The period search procedure is described in detail 
in lScholz fc Eisloffell ll2004ar) and has been thoroughly t ested 
in previous campaigns (IScholz fc Eisloffell l2004bl . l2005h . 

From the 113 VLM members in IC4665, 16 show sig- 
nificant periodic variability in run A and 6 in run C. Only 
two of these objects exhibit a photometric period in both 
campaigns. Final false alarm probabilities for the periodici- 
ties are estimated using the 'bootstrap' approach: For each 
candidate, a set of 10000 random lightcurves was generated 
by shuffling the photometry datapoints, while retaining the 
time sampling. The highest peak in the Scargle periodogram 
was determined for all these test lightcurves, and compared 
with the peak height found for the periodicity in the ac- 
tual lightcurve. A FAP< 0.01% indicates that none of the 
10000 random time series resulted in a Scargle peak higher 
than the one in the observed lightcurves. All our periods 
turn out to be highly significant according to this test. Ta- 
ble 13] lists the derived periods and further results from the 
time series analysis. The period errors AP in the table are 
determined following iHorne fc BaliunasI (|l986l ): for compar- 
ison we also derived AP from the FWHM of the peak in the 
Scargle periodogram which typically gives 1-3% uncertainty. 
The amplitudes in the table correspond to the peak-to-peak 
amplitude of a sinefit to the lightcurve, not the lightcurve 
itself, to avoid being affected by the noise in the photometry. 

An assessment of the completeness of the period search 



Table 3. Results of the period search: objects with significant 
periodic variability in runs A and B. (id: identification number 
from Table [2]|; P: period; AP: period error; A: ampUtude (see 
text); FAP: false alarm probabiHty (see text); N: number of dat- 
apoints used in the time series analysis 



ID P (h) AP (h) A (mag) FAP (%) N 



Run A: 
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0.08 


0.009 


03 


91 


31 


2.99 


0.02 


0.019 


< 0.01 


91 


ou 




U.Ul 


04^ 


^ 01 


yo 
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« 79 
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097 


^ 01 
U.Ul 


oy 


40 


16.5 


0.08 


0.047 


< 0.01 


94 


49 


5.36 


0.08 


0.006 


< 0.01 


90 


55 


7.69 


0.06 


0.007 


< 0.01 


92 


57 


10.3 


0.06 


0.022 


< 0.01 


95 


59 


14.0 


0.04 


0.045 


< 0.01 


93 


63 


9.56 


0.23 


0.015 


< 0.01 


88 


76 


4.26 


0.02 


0.052 


0.05 


92 


88 


3.14 


0.02 


0.012 


0.07 


93 


100 


5.99 


0.02 


0.010 


< 0.01 


93 


110 


9.16 


0.04 


0.033 


< 0.01 


93 


118 


4.49 


0.06 


0.122 


< 0.01 


92 


119 


5.81 


0.11 


0.005 


< 0.01 


84 


Run C: 












28 


4.56 


0.02 


0.011 


0.03 


86 


39 


9.20 


0.03 


0.020 


< 0.01 


73 


57 


16.3 


0.12 


0.022 


< 0.01 


87 


59 


14.4 


0.04 


0.072 


< 0.01 


86 


73 


27.0 


0.35 


0.028 


< 0.01 


87 


75 


4.84 


0.02 


0.033 


< 0.01 


87 



can be made based on the time sampling in our observ- 
ing campaigns. The upper frequency limit for the period 
search, thus the lower period limit, is frequently given as 
I'max = 1/2^, where S is the minimum separation between 
two datapoints (e.g. iRoberts et ai]|l987f ). For our WFI time 
campaigns, this gives a minimum period of Pmin = 0.4 h. 
Since we did not find any period with P < 3 h, we are com- 
plete at short periods. With regular sampling, the longest 
period than can reliably be detected is roughly given by the 
total length of the time series (from first to last datapoint), 
which is 12-13 days in the WFI runs. Nights with very few 
datapoints (sC 5, see Table [1} are unlikely to extend the 
baseline; excluding those nights gives Pmax = 1.3 d for run 
A and 6d for run C. 

For a more realistic estimate of the sensitivity, we 
carried out simu l ations following the recipe given in 
IScholz fc Eisloffell (|2004bl ). In short: Non- variable objects 
were selected from our database. Sine-shaped periodicities 
with periods ranging from 1 to 200 h were added to their 
lightcurves. The noise level of the test objects and the ampli- 
tudes of the artificial periods were chosen to approximately 
reproduce the typical properties of the periods listed in Ta- 
ble[3| We then applied the algorithm described in Sect. l4.l1 to 
these lightcurves and aimed to recover the imposed periods. 

Fig. [5] shows the results fo these simulations for runs A 
and C. The deviation between the imposed period and the 
detected period is negligible for P < 10 h, confirming that 
our period search is optimal for this range of periods. For 
P = 10 . . . 90h there are certain period ranges (e.g., around 
24 h) for which we are unable to obtain reliable detections. In 
many cases, these ranges are narrow (typically 5P < 0.5 h), 
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Figure 3. Phased lightcurves for the periodic objects from run A. Object id and photometric period are indicated. 



P=4 56h Nr. 39 (111 



P=16.3h Nr 59 (111 



particularly for run C. The plots clearly demonstrate that 
there are wide ranges of periods > lOh in which our time 
series analysis is sensitive. For example, in run A we recover 
most periods with P = 40 ... 60 h, and in run C we reli- 
ably recover P = 50 . . . 62 h and P = 65 . . . 90 h. There are 
only negligible windows of sensitivity for P > 95 h, which 
we identify as our upper period detection limit. The non- 
existence of any significant periods with P > 30h in our 
candidate lightcurves thus points to a genuine lack of such 
objects in IC4665 (see Sect.[S}. 



4.2 The TLS lightcurves (run B) 



Nr. 73 [ID ,P=^,'^-Oh Nr. 75 (111 P=A.M^ 



While the photometry from runs A and C comes with simi- 
lar noise characteristics, run B suffers from a comparatively 
: . . : -»»5 - . ■ . ■ - small dynamic range and increased noise level, see Fig. [T] 

- .-. .-. This hampers a direct comparison between the three runs. 

I •■ v ■ . ■■■ v ■ . I .■ ■ . • .■ ■ . Instead of carrying out a time series analysis for all can- 

s' -." ■ ■ I didate lightcurves from run B, we use these time series to 

„„,:•■"'.■ ■ verify the periods reported in Table |3] 

' Out of 20 objects which are periodic in at least one of 

° " pji.. 15 ! a OS ^ i__ IB ! ^jjg \YPI runs, 15 are bright enough to be detected in the 

TLS images. Nine out of these 15 do not show a significant 
Figure 4. Phased lightcurves for the periodic objects from run period according to our tests. For the two objects #55 and 

C. Object id and photometric period are indicated. ^39 confirm the periods listed in Table O within ±0.1 h. 

For two objects we find periods at odds with the ones de- 
tected in the WFI runs: Object #63 has a 7.7 h period (run 
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Period (h) 
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Figure 5. Results of the simulations to test the sensitivity in 
the period search. Plotted is the difference between imposed and 
detected period vs. the period for runs A (upper panel) and C 
(lower panel). The dotted line corresponds to a deviation of 10%. 



A: 9.6h), object #100 a ~ llOh period (ran A: 6.0h). The 
latter case could indicate a long-term trend in the lightcurve 
on timescales of > 4d and thus not detectable in the WFI 
runs. For both ambiguous cases, we re-checked the WFI pe- 
riods and found them to be more convincing. 

Of particular interest are the two objects for which both 
WFI runs yielded a period detection. For candidate #57 
the best periods from the WFI campaigns are 10.3 h in run 
A and 16.3 h in run C. In both cases, however, the peri- 
odograms (Scargle and CLEAN) do exhibit a weaker, but 
still significant peak at the period found in the other run, 
respectively. The two possible periods are neither harmonics 
(Pi = nP2) nor beat periods {1/Pi = I/P2 + n). In run B, 
this object is again variable with a period of 18.4 h, which 
is a clear beat period of 10.3 h, the period from run A (with 
n = —1). Thus, P — 10.3 h is the only period consistent with 
all three time series and is adopted as the most likely rota- 
tion period for this object in the following sections. The fact 
that the most significant period in dataset C differs from our 
consensus period, illustrates the difficulties of assigning pe- 
riods to individual objects based on a single monitoring run. 
Due to the incomplete sampling in ground-based data, the 
lightcurve for a given object can sometimes be fit by more 
than one convincing period. This can be overcome with mul- 
tiple observing runs, as demonstrated here. 



Object #59 shows a period of 14.0 h in run A and 14.4 h 
in run C. The candidate is clearly variable in run B: Its 
lightcurve shows a gradual brightening of ~ 0.04 mag over 
the five nights. In addition, in night 5 the lightcurve ex- 
hibits several outlying datapoints which are 0.05-0.2 mag 
brighter than the average trend, but do not show the typ- 
ical characteristics of a flare (see Sect. |4]). For comparison, 
the average noise level for objects with similar brightness 
is 0.015 mag. However, we do not find any significant pe- 
riod in the lightcurve. In particular, the 'consensus' period 
of ~ 14 h obtained from the WFI lightcurves is not detected. 
Thus, while the variability of this object is confirmed, the 
period is not. 



5 MAGNETIC SPOTS ON VERY LOW MASS 
OBJECTS 

We found periodic variability for 20 likely VLM members 
of IC4665. The best explanation for this type of variability 
is the presence of surface features, co-rotating with the ob- 
jects and thus modulating the flux, as it is usually assumed 
in monitoring studies for this type of targets. Alternative in- 
terpretations include pulsations and eclipsing binaries, but 
both are unlikely for our sample: a) Eclipsing binaries are 
simply too rare to be responsible for a signiflcant fraction of 
our periodic objects, b) There is no evidence which points 
to pulsation in young M dwarfs with the periods observed in 
this study. All further discussions are based on the assump- 
tion that with the perodicities we indirectly observe surface 
features. 

In our observations, the surface features causing the pe- 
riodicities are almost certainly magnetically induced cool 
spots, simply for the lack of an alternative plausible sce- 
nario. It has been suggested that the condensation of dust 
molecules will form 'clouds' in substellar objects, which 
could in principle account for periodic variability as well. 
Our young targets, however, are not yet cool enough for 
this to happen. Models predict t hat dust c ondensation be- 
comes relevant at TeS < 2500 K jHeUing et al. 2008,') . while 
the coolest objects in the sample given in Tableware likely 
to have > 2700 K. Moreover, our targets are expected to 
be mid to late M dwarfs, a spectral range known for strong 
magnetic activity. Thus, by analysing the lightcurves from 
our monitoring campaigns, we can put limits on properties 
of magnetic spots on VLM objects. 

It is useful to recapitulate the conditions required 
to detect a photometric rotation period. Amplitudes and 
lightcurve shape will be influenced by the distribution of the 
spots on the surface. A degree of asymmetry in the distribu- 
tion is required to cause a photometric period, but a slight 
deviation from symmetry would be sufficient to explain the 
small amplitudes seen in this study. Moreover, low-latitude 
spots are more relevant than polar spots for the lightcurve 
appearance. 

The spot coverage (usually measured as filling factor) 
and temperature contrast (difference between photospheric 
and spot temperature) affect the amplitude of the photo- 
metric variability and th us the 'success rate ' in monitoring 
campaigns. As shown in IScholz et all (|2005l ). both param- 
eters are degenerate and cannot be determined separately 
based on single-filter lightcurves. Under the limiting assump- 
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tion that the spot temperature is T = OK, we can define a 
lower limit for the spot coverage: In that case, amplitudes 
of 1-10%, as seen in our period sample, would require a spot 
filling factor of at least 1-10%. In reality, the temperature 
contrast between photosphere and spots is mor e likely to 
be in the range of 100-1000 K (jScholz et al.ll2005l l. implying 
more spot coverage than in the limiting T = OK case. 

Finally, the inclination of the rotation axis against the 
line of sight will affect the lightcurve. For objects seen face- 
on, we expect flat lightcurves, independent of the spot prop- 
erties. Assuming that the orientation of rotation axes in 
open clusters are randomly distributed, this will lower the 
fraction of objects with observed periodicity by a constant 
factor. 

5.1 Spot activity vs. mass 

The periods measured in Sect. 14.11 provide us with two in- 
dicators to assess the level of spot activity in our sample, 
the amplitudes and the fraction of objects with periods. 
The photometric amplitudes found in our period sample are 
^ 0.05 mag in the clear majority of the cases (20 out of 22), 
indicating that VLM objects in IC4665 show only low-level 
photometric variations. This is also reflected in the rms of 
their lightcurves, as discussed in Sect. |3] (see Fig. [T|. As- 
suming that the contamination rate is similar in the total 
candidate sample and in the periodic sample, the fraction 
of objects with periods is 18% (20 out of 113). This num- 
ber increases to 25% (20 out of 80), if we assume that the 
periodic objects have zero contamination (see Sect. [3]). 

These two results will now be compared with the out- 
comes of variability studies for more massive stars to test 
if spot activity changes with object mass. As comparison 
samples we mostly focus on four clusters with similar age, 
IC4665 itself, IC2602, a Per, and the Pleiades. Monitor- 
ing campaigns in these clusters were mostly carried out in 
the V-band, instead of the I-filter we have used. Assum- 
ing typical properties of cool spots, the photometric ampli- 
tudes are likely to be a factor of 1.5-2 larger in the V-band 
(jScholz et al.ll2005l l. This has to be taken into account when 
comparing am plitudes from differ ent authors. 

In IC4665 ]Allain et al.l (j 19961 ) monitored 15 F-K dwarfs 
and find 8 with period with amplitudes 0.01-0.1 mag, a 'de- 
tection rate]_of_53%, significantly higher than in our cam- 
paigns. iBames^eF^l] lfl999l observed all spectroscopic mem- 
bers of IC2602 and measured periods for 29 out of 33 F-K 
dwarfs, a fraction of 88%. With one exc eption, all amp l itudes 
are ^ 0.08 mag. In the a Per cluste r iProsser et al.l (|l993l ) 
and lO'dell fc Collier Cameronl l|l993l ) find photometric pe- 
riods for ~ 90% of the monitored F-K stars with typical 
amplitudes ^0.1 mag. In all these studies, the range of am- 
plitudes is similar to our lightcurves, particularly after con- 
verting from V- to I-band. The detection rate, on the other 
hand, is 50-90% and thus clearly higher than in our VLM 
sample. 

The main observational selection effects that might in- 
fiuence the period detection rates are photometric accuracy, 
sampling of the time series, and a biased choice of targets. 
The photometric accuracy in our campaigns (particularly 
runs A and C) is better than in most aforementioned litera- 
ture studies. In addition, most of the literature studies rely 
on fewer datapoints and are probably less sensitive to short 



periods. Combining both effects, we would expect to find a 
larger fraction of periodic objects, not the opposite. 

All our targets are selected from photometry alone and 
are thus unbiased with respect to the rotation/activity prop- 
erties. To o ur knowledge , the s ame a pplies to the s ample s 
observed in lAUain et all l|l996l ) and ISarnes et al.l l|l999l ). 
The a Per sa mples might be biased towards fast rotating 
objects, since lO'deU fc CoiiteT Camero 3 (|l993l) selected pre- 
sumed fast rotators for their monitoring study. Since the 
rotation- activity relation for pe riods of 1-5 d is fiat in these 
clusters jPatten fc Simonlll996l) . this does not imply a bias 
in activity levels. Overall, selection effects do not appear to 
be responsible for the low number of stars with period in 
our sample. More likely, this finding is caused by a mass 
dependence in the spot properties. 

Monitoring s tudies in NGC2264 ('Lamm' '2003|: 
Herbst et a"l]|2007^ and in the Pleiades (Scholz & Eisl6ff3 
2004bl ) have pointed at a drop in photometric amplitudes 
by a factor of about 2 in the very low mass regime. This has 
been interpreted as evidence for a change in spot properties 
at M ~ 0.3 — 0.5 M0. We cannot confirm this drop in 
amplitudes in IC4665, which may be due to the small 
sample size. Indirectly, however, we confirm these findings, 
as low amplitudes could prevent period detection. 

Possible explanations for low amplitudes and/or a low 
fraction of periodic lightcurves are a) very few spots, b) low 
contrast between spots and photosphere, c) many symmet- 
rically distr ibuted spots o r one polar spot. From multi-filter 
monitoring IScholz et al.l (|2005| ) conclude that a) and /or c) 
are the most plausible scenarios. We note that the current 
results from Zeeman Doppler imaging point towards large- 
scale, poloidal fields for fu lly convective very low mass ob- 
jects (|Morin et al.l l2008bl ). This indicates the presence of 
large polar spots and a lack of spots at low latitudes, which 
reinforces hypothesis c). 

5.2 Long-term spot evolution 

With our three monitoring campaigns covering three years, 
we are for the first time able to probe the long-term evolution 
of magnetic spots in VLM objects. Of particular relevance 
is the comparison between the two WFI time series (run A 
and C): Having similar noise characteristics and sampling, 
they allow us a direct check for long-term variations in the 
lightcurve. 

The most important finding here is that out of 20 ob- 
jects with detected period, only two show a period in both 
campaigns. The lightcurves of those two objects will be fur- 
ther discussed in Sect. 15.31 In the majority of the cases, the 
characteristics of the spots have changed sufficiently after 
three years, preventing the confirmation of the period. Per- 
sistent spot activity is thus a rare phenomenon in young 
VLM objects. 

This result can again be compared with the properties 
of more massive stars. Large samples of solar-mass stars at 
ages 1-5 Myr have been monitored repeatedly in the past, 
allowi ng us to put c onstraints on the stability of the spot ac- 
tivity. ICohen et al.l (|2004 ) have published a five-year study 
of T Tauri stars in IC348. I n agreement with other studies 
(e.g. iNordhagen et al]|2006l ). they find variations in ampli- 
tude and lightcurve shape between different seasons. In their 
sample, 14 out of 27 stars with periodic variability exhibit 
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Table 4. Results from the period fitting to the lightcurves for ob- 
jects #57 and #59. The peri o d erro r AP (h) has been determined 
following iHorne &: BaliunasI jlQ Sd). The period fit did not yield 
significant results for the run B lightcurves; the periods from that 
run are considered unreliable. 



ID 


run 


P (h) 


AP (h) 


A (mag) 


57 


A 


9.99 


0.04 


0.026 




B 


10.45 




0.017 




C 


9.69 


0.04 


0.017 


59 


A 


14.04 


0.04 


0.047 




B 


13.17 




0.027 




C 


14.45 


0.04 


0.083 



the period in all five seasons; 23 of 27 stars in at least two 
out of five seasons. This is consiste nt with results from work 
in the Orion Nebula Cluster, where iRebulil l|200ll ) finds that 
about half of the periods measured in 1996/97 can be recov- 
ered in a dataset obtained about a year earlier, although this 
earlier run had more 'patchy' sampling. Spot stability over 
several year s has also been i nferred from Doppler imaging 
studies (e.g. lVogt et al.lll999l '). 

Thus, for young solar-mass stars about half of the ob- 
jects retain a detectable period over more than one season. 
In contrast, only 10% of the VLM objects in IC4665 exhibit 
a photometric period in both WFI campaigns separated by 
three years. This again points to a change in the spot proper- 
ties in the VLM regime, as discussed in Sect. 15. Il In very low 
mass objects, spot activity causing photometric variability 
on the percent level is rare and transient. Since the ampli- 
tudes in our sample are small and in many cases close to the 
photometric noise limit, only a slight change in lightcurve 
amplitude, caused for example by a re-configuration of the 
spots, could lead to a non-detection in the period search. 
Thus, the lack of objects with stable periods may be a di- 
rect outcome of the drop in amplitudes in the VLM regime. 



5.3 Objects with persistent periodic variability 

The two most interesting objects in our sample are proba- 
bly candidates #57 and #59, both with persistent periodic 
variability in all three observing campaigns, spanning three 
years. With masses ~ 0.2 A/q these two objects belong to 
the lowest mass stars known to date with long-term photo- 
metric variability. They are thus prime targets for follow-up 
studies, e.g. multi-filter monitoring or Doppler imaging. 

We re-determined the periods for #57 and #59 by fit- 
ting sine functions to the lightcurves from each run. The 
free parameters in the fit are phase, amplitude, and period, 
but the period was restricted to a small range around the 
assumed rotation period of 10 and 14 h, respectively (see 
Sect. 14.2) . The fit results are summarised in Table [3] We do 
not attempt to fit all three seasons with a single periodic 
function, because the period uncertainty accumulated over 
a year or more causes a phase uncertainty larger than the 
period itself. In Fig. |6]we show for both objects the three 
lightcurves phased to the period given in Table|4] As pointed 
out earlier, the rotation periods are not recovered in run B, 
but for object #57 the beat period at 18.6 h is detected. The 
period fit for run B are thus not considered trustworthy. 



In the following, we will discuss change in amplitude, 
lightcurve shape, and period over the three seasons. For ob- 
ject #57 the photometric amplitude is constant within the 
errorbars over all three seasons. For object #59, however, 
the amplitude shows striking variations: It drops by a fac- 
tor of about two from 1999 to 2001 and increases then by a 
factor of about three between 2001 and 2002. The two most 
likely scenarios for such changes are a re-configuration of the 
magnetic spots and a strongly variable spot filling factor. 

For both objects, the shape of the lightcurves are in 
good agreement with the sine function. This is not particu- 
larly surprising since the periods have been determined by 
trying to match a sine curve. The more relevant result is 
that the shape of the lightcurve does not change significantly 
over the three seasons. For object #59 this is demonstrated 
in Fig. [7] We plot the datapoints from run A and run C 
in phase to the period measured for run C (14.45 h) scale 
the lightcurve from run A by a factor of 1.766 to match the 
amplitude of run C, and overplot them with the lightcurve 
measured in run C. 

As can be seen in this plot, both sets of datapoints 
cover the same regions in phase space. The deviations be- 
tween both lightcurves are fully consistent with the expected 
scatter of about ±0.015 mag, given as solid lines in Fig. [T] 
This noise estimate includes the combined photometric noise 
of both campaigns (~ 0.01 mag), uncertainty in the period, 
and possible other sources of variabiliy (e.g., microfiares). 
Thus, within the uncertainties the lightcurves from run A 
and C are identical, except for a clear change in amplitude. 
Since any re-configuration of the magnetic spots is expected 
to alter the shape of the lightcurve, this favours a variable 
spot filling factor as explanation for the change in amplitude, 
maintaining the magnetic configuration constant. Assuming 
that spot temperatures are on average constant, changes in 
the spot coverage by a factor of 2-3 on timescales of years 
are required to match the observations. 

Magnetic cycles, as seen in solar-type stars, are one pos- 
sible reason for a long-term variation in the spot filling fac- 
tor. As of today, simulations of the magnetic field generation 
do not find strong evidence for the presence of cycles in fast 
rotating fully connective objects (IChabrier fc Kiikerl l2006l : 
iDobler et al]|2006l ; iBrownind [200a ) . Our result might pose 
a challenge to these models and certainly motivates contin- 
ued monitoring of variable VLM stars. 

In principle, our long-term study also allows us to probe 
differential rotation: If the rotation rate is a function of lat- 
itude, we expect to measure small changes in the rotation 
periods in different seasons, depending on the latitude of the 
dominant spot group. The rate of differential rotation can 
then be estimated simply from the scatter of the measured 
rotation rates Auj. Our candidates #57 and #59 show period 
changes of 0.3 and 0.4 h between run A and C. This trans- 
lates to Auj of 0.4radd~^ for object #57 and O.Sradd^^ for 
object #59, corresponding to 2-3% of the rotational velocity. 

The formal errors from the period fit are smaller than 
that (0.04 h), which would indicate a detection of significant 
differential rotation in those two objects. A more detailed 
look on the lightcurves, however, shows that the actual un- 
certainties in the period measurement are likely to be larger: 
Due to the lack of datapoints outside two consecutive nights, 
the periodogram peaks for run A are broad enough to allow 
for period uncertainties of ~ 0.5 h. More specifically, the pe- 
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Figure 7. Combined datapoints from run A (plusses) and C 
(crosses) for object #59. Both lightcurves are plotted in phase 
to the period measured in run C (14.45 h). The hghtcurve from 
run A has been scaled by a factor of 1.77 to match the amplitude 
in run C. The two solid lines give the range expected from the la 
photometric noise. 



riods obtained for run C give a convincing (albeit not the 
best) fit to the lightcurves from run A. Therefore, we con- 
sider the periods from run A and C to be consistent for both 
objects, and treat the period differences as upper limits for 
the rate of differential rotation. 

From theoretical ground, fully convective objects are 
expe cted to show only little, if any, differential rota- 
tion (|Chabrier fc Kiikeiji2006l: iDobler et ai]|2006l : iBrownind 
|2008| ). Monitoring camp aigns on T Tauri s t ars have con- 
firme d this prediction l|Cohen et al.l l2004l : iHerbst et al.l 
I2OO6I ). with period changes rarely exceeding a few percent. 
For VLM objects there is growing evidence that the differen- 
tial rotation rate is negligible, as found by Doppler imaging 



(< 0.05radd~\ iBarries et"ai1l2005h and Zeeman Doppler 
Imaging f~mradd~^. lMorin et al.ll2008al lbl). Our findings in 
IC4665 are certainly consistent with the expectations. Given 
the lack of objects with persistent variability in this mass 
regime, it is challenging to derive more stringent limits on 
differential rotation from photometric monitoring alone. 



6 ROTATION PERIODS IN IC4665 

As established in Sect. [S] the observed photometric periods 
correspond to the rotation periods of our targetslf] The sam- 
ple of 20 new rotation periods for likely VLM members in 
IC4665 increases the number of objects with known periods 
in this cluster by a factor of three. In Fig. [S] we plot the 
total sample of published periods in IC4665 as a function 
of mass, inc l uding our own periods and the results from 
lAllain et a l.' (1996) for F-K stars in this cluster. We note 
that the AUain et al] (|l996l ) variability campaign is likely 
not highly sensitive to very short periods (< 10 h), due to a 
small number of datapoints per observing night. 

The main feature in Fig. |S] is the cumulation of our 
datapoints at periods < 1 d, or conversely the distinct lack 
of slow rotators in the VLM regime. Our period search may 
be incomplete, but it is definitely sensitive for periods up to 
4d, as discussed in Sect. 14. Il Thus, the observed lack of slow 
rotators in the VLM regime is not induced by our sampling. 
This finding is supported by period samples in other young 
clusters - VLM objects are almost exclusively seen as fast 
rotators, with the exception of objects with ages < 5 Myr, 
when disk-star interactions might still play a role. The data 



^ In rare cases, a particular spot distribution might generate a 
photometric period which is half the rotation period, e.g. when 
two similarly large features are located at opposite latitudes. 
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Figure 8. The total sample of rotation periods in IC4665, plotted 
vs. object mass. Crosses are periods from this pape r for VLM ob - 
jects, triangles show periods from F-K stars ( Allai ~et al.|[l993 ). 
The horizontal line in the top left corner indicates the approxi- 
mate detection limit in our period search. All mass estimates have 
been derived in a consistent way and are thus comparable. 



in IC4665 provide another piece of evidence for a strong 
rotation-mass dependence. 

To date we have no evidence that there is a popula- 
tion of slowly rotating VLM objects not found in the pe- 
riod searches. For the age and mass range considered here, 
photometric amplitudes are observed to be more or less in- 
dependent of rotat ion, without any drop-off at long periods 
ijirwin et al.ll2008bl), which could ca use us to miss long peri- 



ods. We note that 



i llrwin et al.l ||2009D do in fact report such a 
drop-off for stars with 0.5 < M < 1.1 Mq, age 130 Myr, and 
P > 10 d. We do not attribute too much relevance to this 
result in the context of the current study, for three reasons: 
a) The mass range and age are not comparable with our tar- 
gets, b) Our study is focused on the period range < 4 d, for 
which no drop-off in the amplitudes is reported, c) There is 
independent evidence for an upper period limit well below 
10 d in the VLM regime. 

This evidence comes from the available v sin i data ob- 
tained by high-resolution spectroscopy, which is not biased 
with respect to the spot activity. Albeit based on a lim- 
ited sample, the lower envelope of vsmi in the Pleiades 
(120 Myr) for M < 0.4Mb i s weU defined ijTerndrup et all 
I2OOOI ). In lScholz fc Eisloffell l|2004bl ) we have shown that it 
translates into an upper limit for the rotation periods of 
1-2 d, in good agreement with the available periods in the 
Pleiades as well as in IC4665. Thus, the observed lack of 
slow rotators likely is a genuine effect and affects all VLM 
objects, and not only the ones with detected period. 

In the following subsections we discuss rotation in the 
VLM regime. In Sect. I6.II we focus on our new periods and 
evaluate the pre-main sequence rotational evolution based 
on the currently available data. We find clear evidence for 
rotational braking on timescales of 5-50 Myr. In Sect. 16.21 we 
review the currently used wind braking laws for stars and 
brown dwarfs, before applying them to pre-main sequence 
data in Sect. 16.31 Finally, we put these results in context 
and identify problems and future directions for research on 
rotation and magnetic properties (Sect.[7]|. 
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Figure 9. Rotational evolution of VLM objects in the pre- 
main sequence phase. Plotted are periods for objects with 
M < 0.3 Mp) in the clusters a Ori and eOri (average age ~ 
4Mvr. IScholz fc Eisloffeill2004al . lioosi '). NGC2362 (age ~ 5 Myr, 
llrwin et al.l l2008a^ IC466 5 (this work), and NGC2547 (age ~ 
SSMvr. llrwin et al.ll2008d) . The period data for any given cluster 
is plotted at a range of ages scattered randomly around the most 
likely cluster age, for clarity. The large polygons show the median 
period for 0.1 < M < 0.3 Mq. In addition to the groups of clusters 
at 4-5 and ~ 40 Myr, we plot median and 90%-i le for VLM ob- 
jects in NGC2264 at ~ 3 Myr l lLamm et al.ll2005l . M < 0.25 Mq 
on their mass scale). The dashed lines show the period evolution 
assuming angular momentum conservation based on the evolu- 
tionary tracks by iBaraffe et al.l (jl99^ for 0.3 Mq (upper track) 
and 0.2 Mq (lower track). 



6.1 Tracing the pre-main sequence evolution 

There is now a substantial sample of rotation periods for 
VLM objects in young open clusters. With at least 20 
VLM periods in seven pre-main sequence clusters (ONC, 
NGC2264, (J Ori, eOri, NGC2362, IC4665, NGC2547) we 
are in a position to put fir m limits on the rotational brak- 
ing at very young ages (see iHe rbst et al]|2007l :l Irwin et al.l 
In Fig. |9] we plot periods for objects with M < 



0.3 Mq as function of clust er age. As crosses we sh ow pe- 
riods in the clusters a O ri (jScholz fc Eisloffell liooial ). eOri 
(|Scholz fc Eisloffell I2OO5I ) - both plotted at an approximate 
age of 4 Myr -, and IC4665 (this paper, age ~ 36 Myr). 
Periods from the Monitor project are pl otted as small dots , 
for the clusters NGC2362 (age ~ 5 Mvr. llrwin et al.ll2008ah 
and NGC2547 (age ~ 38 Myr, llrwin et al.ll2008bl ). In total, 
we base the following discussion on 72 periods at 4 — 5 Myr 
and 82 periods at 40 Myr. 

For the youngest clusters in this plot, the fraction of ob- 
jects with disks and/or accretion is < 10%, thus rotational 
braking due to disk/star interaction is not relevant here. In 
all five clusters, the object masses have been determined in 
a consistent way, thus the cut-off at 0.3 Mq is consistently 
defined for all ages. For completeness, we also plot the me- 
dian and the 90%-ile of the VL M periods in NGC2264 at 
age of 3 Myr (|Lamm et al.ll2005l ). noting that both masses 
and ages in this cluster are probably not fully comparable 
with the remaining datasets. 

The monitoring campaigns in Orion and IC4665 are 
deeper than the Monitor data, extending into the substel- 
lar regime. Combined with the well-established period-mass 
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relation in the VLM regime, this explains the fact that the 
period samples in these clusters extend to shorter periods. 
The upper period limits, however, are defined by the high- 
est mass objects in those samples and thus comparable. To 
probe the evolution for a typical VLM star, we plot the me- 
dian period for objects with 0.1 < M < 0.3 as large 
polygons, a mass regime where all variability campaigns are 
sensitive. For this purpose we join the samples at similar 
ages, i.e. Orion and NGC2362 (median 37.7 h, based on 57 
objects) as weU as IC4665 and NGC2547 (median 16.5 h, 
based on 78 objects). 

The two dashed lines in Fig. |9] provide the comparison 
with the theory. They show the calculated rotational evo- 
lution assuming angular momentum conservation {Pit) = 
PoX(R{t)/Rof). B oth tracks are based on the radii from the 
iBaraffe et al.l (11998 1) evolutionary tracks for solar-metallicity 
stars. The upper track is for a 0.3 Mq star and can thus be 
compared with the upper period limit, while the lower track 
is for a 0.2 Mq star, meant to provide a comparison with the 
median values. 

Comparing periods and model tracks clearly shows that 
the observations are not in agreement with angular momen- 
tum conservation at timescales of 5-50 Myr. The expected 
median period for a VLM star at ~ 40 Myr is about 8h, 
while the observed median is approximately twice as high. 
Similarly, the upper period limit at 40 Myr is not well- 
represented by the model: The track starts at Po ~ 100 h, 
with only 6 objects (8 ± 3%) have periods exceeding this 
value at young ages. At ~ 40 Myr, however, 24 out of 82 
objects (29 ± 5%) exceed the predicted period of ~ 22 h. 
Excluding the highest 8% of the periods at 40 Myr yields an 
upper limit of about 36 h, a factor of 1.6 too high compared 
with the prediction for angular momentum conservation. 

In summary, rotational braking clearly occurs for ob- 
jects with masses of 0.1-0.3 Mq on timescales of ~ 50 Myr. 
We see an offset between prediction for angular momentum 
conservation and observed periods by a factor of 1.6-2, im- 
plying a reduction in angular momentum by 40-50% com- 
pared with the case without any braking. Note that this 
effect can hardly be explaine d by uncertainties in the radii 
from the iBaraffe et all l| 19981 ) models: We only use the ra- 
tio of radii at different ages, thus any systematic problem 
with the models is likely to be eliminated. Moreover, radii 
from these models are found to be in good agreement with 
observations for e arly /mid M dwarfs (|Mohantv et al..,2003 : 
IScholz et al.1l2007l ). 



In most studies, a value of A'^ = 1.5 is chosen, corre- 
sponding to a mixture between dipole [N ^ 0.5) and ra- 
dial [N = 2.0) field. By adopting a = 1.0 and thus a lin- 
ear relation between B-field and rotation rate, one obtains 
dj/dt oc oj^ and reproduces the Skumanich law u oc t~^^'^ , 
the empirical relationship for the d ecay of rotation r ates of 
G-type stars on the main-sequence l|Skumanichlll973 ). 

This particular choice of the parameter a, however, fails 
to be applicable to the 'ultra-fast rotators' (UFRs) in young 
open clusters. In addition, it does not take into account that 
magnetic activity in dicators are obser ved to saturate at high 
rotation rates (see IPatten fc SimonI [1996, and references 
therein). Therefore, it is usually assumed that the magnetic 
field saturates at a critical rotation rate, i.e. a — 0.0 for 
Lu > cjcrit. With the same wind index as in the Skumanich 
case, Equ. 1 then yields dJ/dt oc o^^rit'^ s-nd J{t) oc exp {t/r), 
with r = l/ti;^^it. By scaling tJcrit appropriately with mass or 
effective temperature, models using this bimodal parameter- 
isation of the angular momentum losses are indeed able to 
reproduce the main features in the availa ble rotation p eriod 
Ce.g. iBouvier et al.l ll997': 'Sills ct al."2000': Irwin ot al. 200"'' 



and V sin i data (e.g. Torndrup ct al. 2000; Rcincrs & Basr 
l2008h . 

Introducing a scaling of the critical rotation rate with 
mass or convective turnover timescale is justified, because 
the magnetic field generation is unlikely to be comparable 
in stars with different interior structure. Main-sequence stars 
with spectral types F to early M possess a radiative core and 
a convective envelope. The convection zone becomes pro- 
gressively deeper with decreasing mass, so that all VLM ob- 
jects are fully convective throughout their evolution. These 
changes are likely to have consequences for the dynamo pro- 
cesses generating the stellar magnetic fields. For example, it 
has been suggested that fully convective objects host an 
dynamo dch abricr & Kiikcr 2006) or a distributed dynamo 
(iDurnev et al]ri993l) instead of the solar-type au) dynamo. 
iBarned (l2003bl ) proposed that the saturated (o = 0.0) and 
linear regimes (a — 1.0) of the angular momentum loss law 
represent objects with 'convective' and 'interface' magnetic 
fields, called objects on the C-sequence or I-sequence. For 
consistency reasons, we will use this nomenclature in the 
following. 



6.2 Wind braking in stars and brown dwarfs 

Stellar winds powered by magnetic fieldsa are the main 
physical process governing the rotational braking on long 
timescales. Angular momentum losses due to stellar winds 
are commonly modeled based on the parameterisation by 
iKawalej l| 19881 ) which is der ived from the analytical wind 
braking law bv iMestell (119841 ): 

Here is the 'wind index' describing the topology of the 
magnetic wind and a is the de pendence of the magn etic field 
on the rotation rate, B oc uj " (|Barnes fc Sofialll996l ). 



6.3 Application to pre-main sequence objects 

Pre-main sequence VLM stars with ages between 10 and 
100 Myr are almost exclusively fast rotators and are thus 
considered to be on the C-sequence with lu > ciJcrit 
(|Pizzolato et al.l |2003| ) . The e-folding timescale for their 
wind braking is proba bly > 1 Gyr (iDelfosse et al.l Il998l : 
IScholz fc Eisloffel |2007| ): the angular momentum losses on 
the pre-main sequence depend primarily on radius with 
to oc R'^^^, according to Equ. 1. This gives an angular mo- 
mentum loss of ~ 30% on timescales of 50 Myr. As eval- 
uated in Sect. 16.11 observations indicate losses of 40-50%. 
Taking into account uncertainties in cluster ages and spin- 
down timescale as well as the fact that a few objects might 
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experience additional braking by star-disk interaction, this 
agrees well with the prediction|j 

On the other hand, the observed estimates of angular 
momentum losses in the VLM regime are nowhere near the 
effect expected for the 1-sequence. Skumanich braking gives 
angular momentum losses of > 70% , i.e. significantly larger 
than the values seen in the VLM regime. This simply con- 
firms that VLM stars on the pre-main sequence have to be 
considered objects on the C-sequence. In summary, the ro- 
tational data for pre-main sequence VLM objects, including 
our new periods in IC4665, fit the current paradigm for the 
rotatio n evolution, as desc ribed in Sect. 16.21 

In lScholz etall (|2007l ). vsini for F to early M stars in 
young stellar associations with ages between 6 and 30 Myr 
have been analysed, a dataset that is complementary to the 
periods shown in Fig. [9] Similarly to the VLM objects, the 
rotational pre-main sequence evolution of solar-mass stars 
is incompatible with Skumanich braking and allows only 
for relatively little angular momentum loss. Thus, in the 
pre-main sequence phase all stars seem to go through a 
C-sequence phase of rapid rotation and weak wind brak- 
ing. They make the transition to the Lsequence and thus 
'switch on' the Skumanich braking at an age tc-i, which 
is a function o f mass. This has also been suggested in a re- 
cent paper bv lMeibom et all (|2009l '). According to the cur- 
rent data, tc -i is 50 - 100 Myr for 1.0 Mq, 0.5-1 Gyr for 
M = O.5M0 (IScholz fc Eislqffelll2007l ). and several Gyrs for 



M < 0.3 M© (|Delfosse et al.lll998l '). The transition between 
I- and C-sequence generates a large gap in rotation rates 
between solar-mass and VLM stars on the main sequence. 



7 TWO REGIMES OF MAGNETIC 
PROPERTIES 

In Fig. [To] we summarise the rotational evolution scheme 
described in the previous Sect. [6] This figure is meant to 
serve as a qualitative illustration and does not provide an 
accurate quantitative description. All stars start on the C- 
sequence with moderate wind braking in the pre-main se- 
quence phase. At a given age tc-i which is a function of 
mass they make the transition to the I-sequence with Sku- 
manich type spindown. In Fig. 1101 we plot the C-sequence in 
solid lines and the I-sequence in dashed lines. For tc-i we 
assume 0.05 Gyr for IMq and 0.5 Gyr for 0.5 M©, consis- 
tent with the currently available constraints (see Sect. 16.3")) . 
Additionally, we assume that the e-folding timescale r on 
the C-sequence scales with object mass. We neglect here the 
dependence of the rotational braking on the radius, which 
is relevant in the pre-main sequence phase. 

With the I- and the C-sequence two regimes of wind 
braking have been identified. In order to explore the under- 
lying physics of these regimes, we look for differences in the 
magnetic properties between C- and I-sequence. The age se- 
quence defined by the transition from C- to I-sequence is 
similarly seen in the fall-off of chromospheric Ha emission, 
which occurs at ^ 50 Myr for G-K stars, at 0.5-1.0 Gyr for 
early M stars (|Hawlev et al.l 0*9991 ) . and at several Gyrs for 



^ Note that the dependence on radius in the wind law has been 
neglect e d in various previous pu blications, e.g. ISchoIz fc Eisloffell 
ll2004lJ ): [He"rbst fc Mundtl ||2005|) . 
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Figure 10. Schematic representation for the rotational evolution 
of stars at three different masses. Solid and dashed lines show 
the spindown on the C-sequence and I-sequence, respectively. The 
vertical dotted lines illustrate the transition from C- to 1-sequence 
for 1 Mq (at 0.05 Gyr) and 0.5 Mq (at 0.5 Gyr). The effect of the 
pre-main sequence contraction has been neglected here (see text). 



Table 5. Qualitative summary of observable properties for ob- 
jects on the C- and 1-sequence. For more references see text. 



C-scq 



1-seq 



rotation 


fast 


slow 


MS wind braking 


oc exp (t)" 


ex y/t 


Hq emission 


yes'' 


no / weak 


rotation/activity relation 


flaf^ 


linear 


differential rotation 


no / weak'' 


strong 


magnetic topology 


poloidal'^ 


poloidal / toroidal 


spot activity 


weak-'' 


strong 



" iBarnesI l|2003b|, and references therein) 
' e.e..| Hawlev et al.l lll 999'. '2000l) 
e.g..|patten fc^ Simo n (1996) _ 



<* Barnes et al. (2005); M orin et al. 



= Donati et al. (■200a ): iMorin et al 
^ this paper 



2008e) 
I 2008al i 



M3-M5 stars (|Hawlev et al.ll2000l ). This may indicate that 
the onset of Hq activity and the transition from C- to I- 
sequence are related, i.e. the 'switching on' of the Skumanich 
braking goes along with a drop in Ha activity. 

As already discussed in Sect. |6]2j the rotation / activity 
relation is flat for very yo ung stars at 1 Myr (Stassu n et al.l 
'2004") and VLM objects (jMohantv fc Basril [20031 ) . both C- 
sequence objects, and linear for slowly rotating solar-mass 
stars on the I-sequence (see also iBarne 3 l2003af l. There are 
clear indications for a change of the magnetic field topology 
at th e transition to fully convective objects (jDonati ct ajj 
I2OO6I ). Additionally, the characteristics of the photometric 
lightcurve, caused by magnetically induced spots, change 
significantly in the VLM regime (see Sect. [5}. Thus, C- and 
I-sequence are likely to constitute two different magnetic 
regimes. We summarise the currently available evidence for 
the bimodality in magnetic/rotational properties in Table [5] 

In the parametrisation for the wind braking given in 
Sect. 16.21 two parameters determine if an object belongs to 
the C- or I-sequence: rotation rate and mass. Since both pa- 
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rameters are interdependent, one of them may be of more 
fundamental relevance than the other. It is particularly 
tempting to explain the whole body of rotation/magnetic 
properties based exclusively on stellar mass, because it de- 
termines the interior structure, as pointed out in Sect. IF!2l 
While VLM objects and very young stars are fully con- 
vective, main-sequence F-K stars harbour a radiative core. 
This difference may be the fundamental reason for the two 

regimes of rotational/magnetic properties^ 

However, as already discussed in IScholz fc Eisloffell 

l|2007l) , the mass and age limits for the transition from C- to 
I-sequence do not coincide with the prediction for the occur- 
rence of the radiative core; thus there is no clear-cut connec- 
tion between the change in internal structure and the change 
in rotation/magnetic properties. For all objects that do 
not remain fully convective throughout their evolution (i.e. 
M > 0.35 Mq), the radiative core in F-K stars is thought 
to appear at < 10-30 Myr and is fully developed after <100- 
200 M yr (jChabrier fc Baraffd Il997l : iD'Antona fc MazziteUil 
1 19941 ). significantly earlier than the transition from C- to 
I-sequence. To maintain the notion that internal structure 
is essential for magnetic properties, one has to assume that 
following the formation of a radiative core, it takes a consid- 
erable timespan before the characteristics of the I-sequence 
(see Table [5]) become observable. 

If interior structure is the underlying reason for the two 
regimes, VLM objects - fully convective throughout their 
evolution - are not supposed to leave the C-sequence. This 
is clearly at odds with the observations. Several studies 
have identified mid/late M dwarfs with very low or unde- 
tected V sin i and acti vity levels below the range expected for 
the s aturated regime (|Hawlev et al.|[2000l : iMohantv fc Basril 
|2003| ). Albeit fully convective, these VLM stars are likely to 
constitute the non-saturated regime, suggesting that slowly 
rotating VLM objects do exist and can operate an I-type 
(Skumanich) rotational braking. From these simple argu- 
ments it seems unlikely that interior structure (or mass) 
alone can account for the two different regimes of mag- 
netic/rotational properties. Thus, two parameters, rotation 
and mass, are required to distinguish between I- and C- 
sequence. 

The fact that we observe two regimes of magnetic prop- 
erties points to a more general problem in the current mod- 
els, as described in Sect. 16.21 The braking laws for I- and 
C-sequence both assume A'' = 1.5, a wind structure some- 
where between dipole and radial field. By holding this pa- 
rameter constant, it is implicitly assumed that all low-mass 
stars and brown dwarfs share the same type of magnetic 
field structure, with saturation on the C-sequence being the 
only difference. As outlined above, there is clear evidence 
for a change of magnetic field properties at the transition 
from C- to I-sequence (see Table [5]), which makes the choice 
of constant questionable. The specific benefit of the value 
N — 1.5 is that is conveniently cancels out the mass loss 
rate in the wind braking law, which is difficult to constrain 
by observations. On the other hand, the explanatory power 
of an angular momentum loss law that does not include the 
mass loss rate is probably limited. 

For a more satisfying understanding of the underlying 
physics, it may be needed to re-investigate the parametrisa- 
tion of the wind braking law in Equ. 1. For solar-mass stars, 
iBarnes fc Sofial (|l996l ') originally proposed two ways to ex- 



plain the ultrafast rotators in young open clusters, objects 
that we would classify as being on the C-sequence: either 
with saturation (a = 0.0 instead of 1.0) or with a low value 
for the wind index until 100 Myr {N < 1, close to a dipole 
field) and then switching to the Skumanich type N — 1.5. 
The second type of parametrisation has not been explored 
in detail and may turn out to be closer to the physical re- 
ality. While saturation of magnetic a ctivity at hi gh rotation 
rates is a well-established fact (see iPatten fc Simon 19961. 
and refere nces therein ), it is probably the filling factor that 
saturates (|Saai|[l99^ ). When choosing a — 0.0 in the wind 
braking law, however, we assume that the magnetic field 
strength is saturated. It may be necessary to resort to mod- 
els with a and N as free parameters, thus accomodating both 
for saturation and a change in the magnetic wind structure. 

The choice of A = 1.5 may not be appropriate for the 
Skumanich type I-sequence either : Recent numerical sim- 
ulations by iMatt fc Pudrit j l|2008l ) indicate a significantly 
smaller value for this parameter (about 0.7 instead of 1.5). 
This would imply that angular momentum losses depend at 
least weakly on the mass loss rate {dj/dt oc M°'^), which 
is probably a more plausible description than one without 
dependence on M. 

From these considerations is becomes clear that the cur- 
rently available theoretical framework for the stellar spin- 
down is lacking a satisfactory description of the physics of 
the rotation-wind connection. Future work should aim to 
verify the bimodal wind braking law as given in Equ. (1) and 
explore possible alternatives. Observationally, the emerging 
problems encourage to obtain more reliable constraints on 
mass loss rates, as this may turn out to be crucial to improve 
the modeling of rotational evolution. Having said this, the 
bimodal schematic used in current models and described in 
detail above has been proven to have considerable predic- 
tive power. Given a proper calibration, it allows us to use 
rotation rates as indirect indiator of stellar ages and masses. 



8 CONCLUSIONS 

We have monitored a large sample of very low mass objects 
and brown dwarfs in the young cluster IC4665 over three 
seasons, aiming to measure rotation periods and to assess 
the long-term evolution of spot activity. For 20 objects pho- 
tometric periods have been derived, which likely represent 
the rotation periods. Two objects show persistent variability. 
For these two, we find a consistent period that is detected 
in all three datasets. However, this is not necessarily the 
most significant period in each dataset. This illustrates that 
periods for individual objects obtained from a single mon- 
itoring campaign might not always represent the rotation 
period. The periods range from 3 to 30 h, with photomet- 
ric amplitudes from a few mmag to 0.12 mag. We find that 
we are sensitive to periods longer than 30 h, i.e. the lack of 
longer periods is not a bias in our dataset. This indicates 
a lack of slow rotators among VLM objects in this cluster. 
We confirm that flares are very rare events in very low mass 
objects. 

Compared with solar-mass stars, the very low mass ob- 
jects in IC4665 show weak photometric variability with low 
amplitudes, which could be due to a change of spot prop- 
erties at the fully convective boundary, for example to a 
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more symmetric spot distribution. Persistent variability is 
rare, again in contrast to the results for more massive stars. 
We do not see evidence for differential rotation. One object 
with persistent periodic variability changes its amplitude by 
a factor of two over the course of 3 years, while the shape of 
the lightcurve remains the same, which could be explained 
by a long-term magnetic cycle. Our sample of lightcurves 
in IC4665 provides a fundament for future studies of the 
long-term evolution of activity in the VLM regime. 

We find that our periods fit into a scenario where 
VLM objects experience moderate wind braking in the pre- 
main sequence phase. They are inconsistent with the strong 
Skumanich-type wind braking. This is in line with the cur- 
rent paradigm for the angular momentum evolution of low- 
mass stars, which includes a bimodal nature of the wind 
braking. We point out that this bimodality is also seen in 
various properties associated with the magnetic activity, in- 
dicating a more fundamental physical difference. While the 
current parameterised description of wind braking is success- 
ful in reproducing main features in the period distributions 
in open clusters, it does not provide a satisfactory physical 
picture. 



ACKNOWLEDGMENTS 

We would like to thank Ansgar Reiners, Sean Matt, 
Jonathan Irwin, and Jerome Bouvier for instructive dis- 
cussions regarding topics related to this paper. This work 
was partially supported by Deutsche Forschungsgemein- 
schaft (DFG) grant Ei 409/11-1 and 11-2. 



REFERENCES 

AUain S., Bouvier J., Prosser C, Marschall L. A., Laakso- 

ncn B. D., 1996, A&A, 305, 498 
Bailer- Jones C. A. L., Mundt R., 2001, A&A, 367, 218 
Barafi'c I., Chabrier G., Allard F., Hauschildt P. H., 1998, 

A&A, 337, 403 
Barnes J. R., Cameron A. C, Donati J.-F., James D. J., 

Marsden S. C, Petit P., 2005, MNRAS, 357, LI 
Barnes S., Sofia S., 1996, ApJ, 462, 746 
Barnes S. A., 2003a, ApJ, 586, L145 
— , 2003b, ApJ, 586, 464 
— , 2007, ApJ, 669, 1167 

Barnes S. A., Sofia S., Prosser C. F., Stauffer J. R., 1999, 
ApJ, 516, 263 

Berger E., 2006, ApJ, 648, 629 

Bouvier J., Forestini M., AUain S., 1997, A&A, 326, 1023 
Browning M. K., 2008, ApJ, 676, 1262 
Chabrier G, Baraffe I., 1997, A&A, 327, 1039 
Chabrier G, Kiiker M., 2006, A&A, 446, 1027 
Cohen R. E., Herbst W., Williams E. C, 2004, A J, 127, 
1602 

DAntona F., Mazzitelh I., 1994, ApJS, 90, 467 
Delfosse X., Forveille T., Perrier C, Mayor M., 1998, A&A, 
331, 581 

Dobler W., Stix M., Brandenburg A., 2006, ApJ, 638, 336 
Donati J.-F., Forveille T., Cameron A. C, Barnes J. R., 

Delfosse X., Jardine M. M., Valenti J. A., 2006, Science, 

311, 633 



Donati J.-F., Morin J., Petit P., Delfosse X., Forveille T., 
Auricrc M., Cabanac R., Dintrans B., Fares R., Gastine 
T., Jardine M. M., Lignieres F., Paletou F., Velez J. C. R., 
Theado S., 2008, MNRAS, 390, 545 

Durney B. R., De Young D. S., Roxburgh I. W., 1993, 
Sol. Phys., 145, 207 

Gelino C. R., Marley M. S., Holtzman J. A., Ackerman 
A. S., Lodders K., 2002, ApJ, 577, 433 

Hallinan G, Antonova A., Doyle J. G, Bourke S., Brisken 
W. F., Golden A., 2006, ApJ, 653, 690 

Hawley S., Reid I. N., Gizis J., 2000, in Astronomical Soci- 
ety of the Pacific Conference Series, Vol. 212, From Giant 
Planets to Cool Stars, Griffith C. A., Marley M. S., eds., 
pp. 252— h 

Hawley S. L., Reid I. N., Gizis J. E., Byrne P. B., 1999, in 
ASP Conf. Ser. 158: Solar and Stellar Activity: Similar- 
ities and Difi^erences, Butler C. J., Doyle J. G, eds., pp. 

63-+ 

Helling C, Woitko P., Thi W.-F., 2008, A&A, 485, 547 
Herbst W., Dhital S., Francis A., Lin L., Tresser N., 

Williams E., 2006, PASP, 118, 828 
Herbst W., EislofTel J., Mundt R., Scholz A., 2007, Proto- 

stars and Planets V, 297 
Herbst W., Mundt R., 2005, ApJ, 633, 967 
Home J. H., Baliunas S. L., 1986, ApJ, 302, 757 
Irwin J., Aigrain S., Bouvier J., Hebb L., Hodgkin S., Irwin 

M., Moraux E., 2009, MNRAS, 392, 1456 
Irwin J., Aigrain S., Hodgkin S., Irwin M., Bouvier J., 
Clarke C, Hebb L., Moraux E., 2006, MNRAS, 370, 954 
Irwin J., Hodgkin S., Aigrain S., Bouvier J., Hebb L., Irwin 

M., Moraux E., 2008a, MNRAS, 384, 675 
Irwin J., Hodgkin S., Aigrain S., Bouvier J., Hebb L., 

Moraux E., 2008b, MNRAS, 383, 1588 
Kawaler S. D., 1988, ApJ, 333, 236 
Kiraga M., Stepien K., 2007, Acta Astronomica, 57, 149 
Lamm M. H., 2003, PhD thesis, PhD Thesis, Com- 
bined Faculties for the Natural Sciences and for Math- 
ematics of the Ruperto-Carola University of Heidelberg, 
Germany (Naturwissenschaftlich-Mathematische Gesamt- 
fakultat der Universitat Heidelberg, Germany). XII -|- 143 
pp. (2003) 

Lamm M. H., Mundt R., Bailer- Jones C. A. L., Herbst W., 

2005, A&A, 430, 1005 
Landolt A. U., 1992, AJ, 104, 340 
Matt S., Pudritz R. E., 2008, ApJ, 678, 1109 
Mcibom S., Mathieu R. D., Stassun K. G, 2009, ApJ, 695, 

679 

Mermilliod J. C, 1981a, A&AS, 44, 467 

— , 1981b, A&A, 97, 235 

Mcstcl L., 1984, in Lecture Notes in Physics, Berlin 
Springer Vcrlag, Vol. 193, Cool Stars, Stellar Systems, and 
the Sun, Baliunas S. L., Hartmann L., eds., pp. 49 — h 

Mohanty S., Basri G, 2003, ApJ, 583, 451 

Mohanty S., Jayawardhana R., Basri G, 2004, ApJ, 609, 
885 

Morin J., Donati J.-F., Forveille T., Delfosse X., Dobler 
W., Petit P., Jardine M. M., Cameron A. C, Albert L., 
Manset N., Dintrans B., Chabrier G, Valenti J. A., 2008a, 
MNRAS, 384, 77 

Morin J., Donati J.-F., Petit P., Delfosse X., Forveille T., 
Albert L., Auriere M., Cabanac R., Dintrans B., Fares 
R., Gastine T., Jardine M. M., Lignieres F., Paletou F., 



16 Scholz, Eisldffel & Mundt 



Ramirez Velez J. C, Theado S., 2008b, MNRAS, 390, 567 
Nordhagen S., Herbst W., Rhode K. L., Williams E. C, 

2006, AJ, 132, 1555 
O'dell M. A., Collier Cameron A., 1993, MNRAS, 262, 521 
Patten B. M., Simon T., 1996, ApJS, 106, 489 
Pizzolato N., Maggio A., Micela G., Sciortino S., Ventura 

R, 2003, A&A, 397, 147 
Prosser C. F., Shetrone M. D., Marilli E., Catalano S., 

Williams S. D., Backman D. E., Laaksonen B. D., Adige 

v., Marschall L. A., Stauffcr J. R., 1993, PASP, 105, 1407 
RebuU L. M., 2001, AJ, 121, 1676 
Reid I. N., Gizis J. E., 1997, AJ, 113, 2246 
Reinors A., Basri G., 2007, ApJ, 656, 1121 
— , 2008, ApJ, 684, 1390 

Roberts D. H., Lohar J., Drolier J. W., 1987, AJ, 93, 968 
Robin A., Creze M., 1986, A&A, 157, 71 
Robin A. C, Reyle C, Derriere S., Picaud S., 2003, A&A, 
409, 523 

Rodn'guez-Lodcsma M. V., Mundt R., Eisloffcl J., Herbst 
W., 2009, in American Institute of Physics Conference Se- 
ries, Vol. 1094, American Institute of Physics Conference 
Series, Stempels E., ed., pp. 118-123 

Saar S. H., 1996, in lAU Symposium, Vol. 176, Stellar Sur- 
face Structure, Strassmeier K. G., Linsky J. L., eds., pp. 
237-+ 

Scargle J. D., 1982, ApJ, 263, 835 

Scholz A., Coffey J., Brandeker A., Jayawardhana R., 2007, 

ApJ, 662, 1254 
Scholz A., Eisloffel J., 2004a, A&A, 419, 249 
— , 2004b, A&A, 421, 259 
— , 2005, A&A, 429, 1007 
— , 2007, MNRAS, 381, 1638 

Scholz A., Eisloffcl J., Froebrich D., 2005, A&A, 438, 675 
Sills A., Pinsonneault M. H., Terndrup D. M., 2000, ApJ, 
534, 335 

Skumanich A., 1972, ApJ, 171, 565 
Spiegel E. A., Weiss N. O., 1980, Nature, 287, 616 
Spruit H. C, van Ballcgooijcn A. A., 1982, A&A, 106, 58 
Stassun K. G., Ardila D. R., Barsony M., Basri G., Mathieu 

R. D., 2004, AJ, 127, 3537 
Stetson P. B., 1987, PASP, 99, 191 

Terndrup D. M., Stauffer J. R., Pinsonneault M. H., Sills 
A., Yuan Y., Jones B. F., Fischer D., Krishnamurthi A., 
2000, AJ, 119, 1303 

Vogt S. S., Hatzes A. P., Misch A. A., Kiirster M., 1999, 
ApJS, 121, 547 



